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The potentials and fields in a two-dimensional model of a charge- 
coupled device (CCD) are studied. We assume no mobile minority car- 
riers have been injected into the CCD and that the electrode voltages do 
not vary with time. The nonlinear equations describing the devices are 
first linearized using the depletion layer approximation. The linearized 
equations are then solved approximately by a fitting technique. Both 
surface and buried channel CCD's are considered. The accuracy and cost 
of obtaining the solution is discussed. This work is a continuation of a 
study initiated in an earlier paper. 1 

I. INTRODUCTION AND SUMMARY 

In this paper we study the electrostatic potential distribution and 
fields in a two-dimensional model of a charge-coupled device 23 (CCD). 
This work is a continuation of a study initiated in an earlier paper 1 
hereafter referred to as I. In I we considered a static, two-dimensional 
model with no mobile charge, and with electrodes so close together 
that they could be assumed to touch. We showed there that the deple- 
tion layer approximation 4 could be used to linearize the potential 
equations, and the linearized equations were then solved analytically. 
The numerical evaluation of these solutions was shown to be very 
accurate and cheap. 

We extend the model of I to allow for gaps between the plates. 
Our purpose here is twofold. We want to examine the dependence 
of the potentials and fields in a CCD on various design parameters. 
As we show, our model allows considerable flexibility in describing 
various electrode configurations. In addition, however, we want to 
demonstrate a method of numerically solving the potential equations 
which we believe is of considerable interest in itself. 

Both surface 2 and buried channel 56 CCD's are considered. However, 
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as in I, only the analysis for buried channel CCD's is given. The results 
for a surface CCD can be obtained as special cases of the results for 
buried channel CCD's. We refer the reader to I for a more detailed 
derivation of the equations and for a discussion of their linearization 
by the depiction layer approximation. 

As we show by examples, fairly complicated models of CCD's can 
be analyzed at moderate cost by the methods of this paper. Neverthe- 
less, the cost of using the methods of I to analyze a CCD with zero 
separation between the electrodes is typically an order of magnitude 
less than the cost of using the methods of this paper to analyze a CCD 
with nonzero electrode separation. This suggests that, in any com- 
plicated design problem, the methods of I should be used to rough out 
a solution, and then the solution should be "fine tuned" by using the 
methods of this paper. In addition, it is shown that, when the gaps 
between the electrodes are of the order of 1 jum, the potentials of 
interest are approximated well by the potentials in the same CCD 
with the electrode separation set to zero. 

The nonlinear equations and boundary conditions defining the 
boundary value problem are introduced in Section II. The linearized 
equations arc also introduced there. In Section III we discuss in some 
detail how we obtain approximate solutions to the linearized problem. 
The reader uninterested in the mathematical details should skip Section 
III and proceed directly to Section IV, which is devoted to examining 
some of the solutions with emphasis on how they are affected by 
changes in the design parameters. We examine a number of different 
design parameters, particularly for buried channel devices. The accu- 
racy of the solutions and the cost of obtaining them is considered in 
Section V. Finally, some mathematical details are contained in two 
appendixes. 

II. THE POTENTIAL EQUATIONS 

We consider CCD's in which the minority carriers are holes and the 
underlying substrate is n-type silicon. The analysis can be modified in 
an obvious way to describe the case where the minority carriers are 
electrons and the substrate is p-type silicon. 

A buried channel CCD consists of a substrate of n-type silicon on 
top of which there is a layer of p-type silicon. The p-type layer is 
covered with a layer of Si0 2 , and closely spaced electrodes are placed 
on top of the oxide layer. A schematic diagram of such a device is 
shown in Fig. 1 with some typical dimensions indicated. A surface CCD 
is the same, except that the p-type layer is missing. 
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Fig. 1 — A schematic diagram of a buried channel CCD. 



We study here the static potential and fields in either a surface or 
buried channel CCD in the absence of mobile charge. Since the length 
in the z-direction of each plate is much greater than its width in the 
z-direction, near the center of the plates (z=0) the field is essentially 
two-dimensional; therefore, we treat the problem as two-dimensional. 

We assume the bottom (n-type) substrate is infinitely thick. The 
field can penetrate into the substrate little beyond a depletion depth 
and, since for typical voltages the depletion depth ranges from 7 to 
20 /im and the thickness of a typical device is 100 ^m, this is a very 
good approximation. 

It is assumed that there are gaps between the electrodes. We also 
make the approximation that the electrodes have zero thickness. 
Although this is a rather drastic simplification, we feel the essential 
effects of the gaps between the electrodes are still properly described. 
It will be seen later that electrodes of rectangular cross section could 
be studied, although at much greater cost. We further assume that the 
medium surrounding the electrodes has the same dielectric constant as 
the Si0 2 . This is very reasonable, since in practice a CCD is covered 
with a dielectric coating. Two basic types of metalization are studied, 
single level and double level. In double-level metalization, two layers 
of electrodes are separated by an oxide layer. This is illustrated 
schematically in Fig. 2. We simulate this situation by assuming that 
the potential distribution in the gaps between the electrodes in the 
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Fig. 2— A schematic diagram of one cell of a buried channel CCD with double- 
level metalization. 

upper layer of electrodes is a known function. (Typically, the potential 
is assumed to vary piece wise linearly.) In single-level metalization, the 
upper level of electrodes is missing. Here we assume the dielectric 
coating over the electrodes is infinitely thick. Again, this is a reasonable 
assumption, since typically the field will have died out before reaching 
the surface of the dielectric coating. 

Finally, we assume the structure to be periodic in the z-direction, 
which in the usual mode of operation is an excellent approximation. 

The boundary value problem corresponding to our model of a buried 
channel CCD can be described by a system of partial differential 
equations which we wish to write in terms of dimensionless quantities. 
All dimensional quantities (measured in rationalized MKS units) will 
be starred, with the exception of a few obvious physical parameters. 
Corresponding unstarred quantities will be dimensionless. The physical 
parameters of the problem are ei and e 2 , the permittivity of the oxide 
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and silicon, respectively; —e, the charge of an electron; Boltzman's 
constant k; the absolute temperature T; the length of a unit cell in the 
device L*; the separation of the two layers of metalization h-u and 
h\ and h 2 , the thickness of the oxide layer and the p-layer, respectively. 
The donor number density of the n-type substrate is N' D , which in the 
usual method of fabricating a CCD is a constant. However, we assume 
the acceptor number density in the p-layer is given by the expression 7 

«<w-c:«p(-(£^jl)«.g)-«. a) 

where C) is the number density of acceptor ions at the upper surface 
of the Si. 

Now define the (dimensional) Debye length X D , 

X = (e 2 kT/e*N' D )K (2) 

Then the dimensionless lengths are defined as 

x = x*/X D) V = V*/\d, L = L*/X D , h a = K/X D , (a = ±1, 2). (3) 

The dimensionless potential is related to the dimensional potential by 

?(s, y) = e<p*(x*, y*)/kT. (4) 

If we set 

C. = C\/N* D , (5) 

then the dimensionless p-layer acceptor density, <r{y), is 



<V) 



= C.ex P |-(L^C. 



-1. (6) 



In the strip ^ x 5| L, let «po denote the electrostatic potential 
above the oxide layer, — oo < y ^ in the case of single-level metaliza- 
tion and — h-i ^ y ^ in the case of double-layer metalization. 
Further, let <pi denote the potential in the oxide layer, ^ y ^ hi', tp% 
the potential in the p-type layer, hi ^ y ^ h 2 ; and <p 3 the potential in 
the n-type substrate (see Fig. 2). Then in the dimensionless form, the 
potential equations are 



VVo = 0, 




y^o, 


(7) 


V^i = o, 




O^y^hi, 


(8) 


V*<pl = a(y) 




hi ^ y ^ h 2 , 


(9) 


V 2 ip 3 = exp 


M -i, 


h 2 ^ y <<*>, 


(10) 



where V 2 is the two-dimensional Laplace operator. The standard 
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electromagnetic boundary conditions are as follows: 

| <po(x, — oo ) | < oo (single-level metalization), and (11) 

ipo(x, —h-i) = Uix) (double-level metalization), (12) 

where U(x) is a given periodic function of period L, assuming on each 
electrode of the second level of metalization the constant voltage of 
the electrode and a specified potential between the electrodes. In 
reality, of course, the true potential in the gaps of the second level of 
metalization is unknown a priori. However, as indicated in Fig. 2, 
typically the semiconductor cannot "see" these gaps since they are 
shielded by the electrodes of the first level of metalization. Thus we 
simulate the exact boundary conditions in the gaps, most often by 
assuming the potential varies linearly from one electrode to another. 
We feel this is a good approximation, since we have performed calcu- 
lations of the potential in the semiconductor with several different 
assumptions about the variation of the potential in the gaps, and the 
results were essentially identical. Further, 

<p (x, 0) = V; = vxix, 0), (j = 1, 2, • • • , p), (13) 

M*, 0) = <Pt(x, 0), ^ (*, 0) = *£ (x, 0) + Pg (x), (14) 

where Vj is the constant voltage of the ith electrode in the first level 
of metalization, eq. (13) holds on each of the p electrodes, and (14) 
holds in the gaps between the electrodes. Typically, p g (x) = 0, but in 
some cases it may describe a deliberately implanted surface charge in 
the gaps. In any event, p g (x) is a known function of x in the gaps, and 
p (x) m on the electrodes. Finally, 

<Pi(x, hi) = <p*(x, h), v ^ (x, h) = ^ (x, ht) + Q(x), (15) 

»(* h 2 ) = *>,(*, h 2 ), ^ (x, h 2 ) = ^ (x, h), (16) 

pifc oo) «=0, (17) 
and for all y 

*(0, y) = *(L, y), % (0, y) = ^ (L, y). (18) 

In (15), Q(x) is a known, periodic surface charge density, which may 
include deliberately implanted charges, 8 and 

v = eiAa- (19) 
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All the boundary conditions (11) to (12) and (14) to (17) hold for 
g x g L. 

The equations for the potential in a surface CCD are essentially the 
same, except the p-type layer is eliminated. We only give the analysis 
for the buried channel CCD. The results for the surface CCD can be 
obtained from those for the buried channel CCD by setting <r(y) = 0, 
hi = h 2 , and <p 2 = <p3. In either case, the fields are obtained from the 
potential by 

E =- Vv>. (20) 

The results of I show that the system of eqs. (7) to (18) can be 
accurately solved by the method of finite differences only at great 
expense for even the simplest of devices. However, it was shown in I 
that, for the simpler problem studied there, the nonlinear boundary 
value problem could be replaced by a linear boundary value problem. 
This linear problem was solved analytically. It was then shown that 
under appropriate conditions the solution of the linear problem was an 
excellent approximation to the solution of the nonlinear problem in 
the p-type layer for a buried channel CCD and near the oxide-semi- 
conductor interface for a surface CCD. The condition for the approxi- 
mation to be a good one is, basically, that the potential along the line 
y = h 2 be large and negative. That condition holds in this problem, as 
we show later by example. Although the linear problem is much more 
complicated in this case because of the gaps between the electrodes and 
we have been unable to solve it analytically, we have been able to 
obtain good approximate solutions of it. For these reasons we now 
formulate the linear problem. 

The linear equations are based on the depletion layer approximation, 
and we refer the reader to I and Ref. 4 for a detailed discussion. The 
linearization consists of replacing the single region h 2 ^ y < <» by two 
regions, h 2 ^ y ^ hi = h 2 + & (the depletion layer) and h 3 ^ y < °° , 
and replacing the single nonlinear equation (10) by a different linear 
equation in each of these subregions (see Fig. 3). 

VVofo V) - 0, y ^ 0, (21) 

VVi(*, V) - 0, £ y £ K (22) 

VV 2 (x, y) = a(y), hi*y£ h 2 , (23) 

VVsfa y) - -1, fa ^ y ^ fa = h 2 + R, (24) 

vvifo y) = Mx, y), fa ^ y < » . (25) 

In addition to \J/ , ip h \f/ 2) and \J/ 3 satisfying boundary conditions (11) to 
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Fig. 3 — A schematic diagram of the depletion layer approximation for one cell of a 
buried channel CCD with double-level metalization. 

(16), we have the boundary conditions for ^ x ^ L 



d^ 



dfa 



M*, h) = Ux, h), ^ (x, h 3 ) = -^ (x, h 3 ), (26) 

Mx, ») =0, (27) 

and the ^ a , (« = 0, 1, 2, 3, 4) all satisfy (18). The pseudodepletion 
depth A is best given by 

£ = _/i + fc 2 _ft 1+ ^+ jY^-Jn+^y -1-27 

' ? T Q " + 2 C 1 (s-^ + v)^)^] 1, (28) 

In (28), * "I 

Q„ =1 j o Q(x)dx 
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and V is the average electrode voltage. As shown in I and Ref. 4, \J/ , 
\f/i, and \f/ 2 are quite insensitive to the choice of R, and an optimal 
choice of R tends to minimize sup [ 4n(x, fh) + 1 1 . 

0<x<L 



III. SOLUTION OF THE LINEARIZED POTENTIAL EQUATIONS 

The method we shall use to solve the linear system of elliptic eqs. 
(21) to (27) has much in common with previous work 9 on the classical 
problem of a single linear elliptic boundary value problem on a simply 
connected domain in the plane. The technique used in Ref. 9 is quite 
simple: Construct a family of particular solutions of the partial 
differential equation and, using a finite linear combination of these 
particular solutions, obtain a Chebyshev fit to the boundary conditions 
at a finite number of points on the boundary of the domain. It was 
shown that, as more and more particular solutions are taken in the 
linear combination and more and more points are chosen in the fit on 
the boundary, the linear combination converges to the true solution. 

In this paper we construct a family of particular solutions of (21) 
to (27) . These solutions depend linearly on a finite number N of param- 
eters and satisfy all the boundary and interface conditions except that 
they do not assume the correct voltages on the electrodes at y = 0. 
We then complete the analogy with Ref. 9 by picking M points on the 
plates, x { , 1 ^ i ^ M , where M ^ N, and force the potential to take 
on the correct value at the points x t in the least-squares sense. 

We obtain the family of particular solutions in the form of Fourier 
series. We assume as given the Fourier series expansions of U(x) and 
Q(x): 

U(x) = *ao + £ G n (x), (29) 

71 = 1 

Qix) = |f + t *.(*), (30) 



where 



and 



G n (x) = a„ cos \„x + /3„ sin \„x, (31) 

$n(x) = f „ cos X n x + £„ sin \ n x, (32) 

X n = (2nr)/L. (33) 



Further, from I, we can write down formal expressions for \p\, \f/ 2 , ^3, 
and \f/i which satisfy (22) to (25) and boundary conditions (15), (16), 
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(18), (26), and (27). Let 

F n {x) = a n cos \ n x + b n sin \ n x, (34) 

E ± = l± l/v, (35) 

A* - 1 =fc (1 + K 2 )*, (36) 

M n (y) = {E+A+ + E-Ane- n ^- h ^}e-^y 

+ {E-Ate- n " hl + £ + A„-e- 2X «^}e x «'', (37) 

L.fo) = 2{A+erX-» + a^x. (*,-,,)}, (38) 

where 77 is given by (19), \ n by (33), and o , a„, b„ (n = 1, 2, • • •) are 
unknown constants. Then these expressions are 

Mx, V) = (A** + B ) + ( c <*o + D)y 

a. f \ F (x) Mn(y) + * (x) Ln(hl) Sinh Ky (39) 

fate, y) = 2 ^°( 1 + a 3 - a 2 ) - (A 3 - a 2 ) 2 

-(00 + 2A 2 - 2A 3 )(?/ - A 2 ) + 2 f" (2/ - £)*($)#] 

+ Z\F n (x) + ^(x)-^-\ w - m , (40) 



Ms, 2/) = i[6o(l + Ai - y) - (y - M 2 ] 



+ E 



F„(aO + *»(*) 



n=l 



sinh X„/ii 

T?Xn 

sinh \ n hi 



L n {y) 
M n (0) 



, (41) 



V^n 



exp [— Vl + \l(y - A 3 ) - X n A 3 ] ti #, 

MJOT ' c } 



where 



A = = ( 1 + A, - A x + 



Ai\ 



(43) 



- i [(A 8 - A 2 )(A, + A 2 - 2A0 + - (fo + 2 A 3 - 2A 2 ) J , (44) 

C=- 1/(27,), (45) 

D = [f + 2(A 3 - A 2 ) - 2 £' *(*)#] /(2u), (46) 
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and 

a = Qa - B)/A. (47) 

It should be noted that 

*i(*, 0) = ha + t F n (x). (48) 

n=l 

In the case of single-level metalization we can write down an expres- 
sion for a solution of (21) in y ^ which satisfies boundary conditions 
(11): 

Mx) = loo + L F n (x)e*»*. (49) 

In the case of double-level metalization, a solution of (21) in 
— h-i^y^0 satisfying boundary condition (12) is 

if \ 1 /, . 2/ \ ,^cw\ sinh \ n (y + ^_0 



-?(£)-£*«= 



sinh X„s, 



sinh X„/i_i 



Note that these solutions have been constructed so that from (48) 
and either (49) or (50), for ^ x ^ L, 

*o(x, 0) = ^(a, 0). (51) 

[Note also that, term by term, (49) is the limit as /i_i — ><*> of (50).] 
Equations (34) to (50) contain expressions for \f/ a , ^ a ^ 4, which 
satisfy the differential equations (21)-(25) and all the required bound- 
ary and interface conditions except condition (13) on the plates and 
the normal derivative condition of (14) in the gaps. These particular 
solutions contain the unknown parameters a , a„, b n , (n = 1, 2, ■ • •), 
which remain to be determined. At this point it might be assumed that 
the series should be truncated at some n = N and the 2JV + 1 coeffi- 
cients a , a n> b n , n = 1,2, • • ■ , .V, be determined directly by making a 
least-squares fit to the remaining boundary conditions. However, it 
can be shown 10 that if x$-\ < Xj are the end points of an electrode on 
y = then for x,-_i < x < x jt and x near x,, say, d\J/ /dy(x, 0) will 
behave like (xj — x)~* plus a power series in (xy — x)K This implies 
that the Fourier series for \f/o(x, 0) converges very slowly. In fact, we 
have found that it is often necessary to take up to 2000 terms in the 
series to represent ipo(x, 0) adequately. This makes it impractical to 
use the Fourier coefficients themselves as the parameters to be deter- 
mined directly by a least-squares process. 
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Instead, we used the following technique. We approximated the 
charge density on y = by a finite sum of known functions 

p(x) = Pa (x) + £ PiPi (x). (52) 

The functions pj(x) are zero except on the electrodes and are of several 
types, as shown in Fig. 4. The function p (x) is a periodic, triangular 
spline for the ends of the device, as shown in Fig. 4a. Corresponding to 
the edge of each plate, there is a discontinuous triangular spline, Figs. 
4b and 4c, and singular splines of the form | x — x k \ - *, Figs. 4d and 4e. 
The remaining Pj{x), whose supports lie wholly interior to the elec- 
trodes, are triangular splines as shown in Fig. 4f. 

Now each unknown parameter a , a n , b n (n = 1, 2, • • •) is deter- 
mined as a linear sum of the N parameters p h 1 ^ j ^ N, by equating 
p{x), given in (52) with d^ /dy(x, 0) - dfa/dyix, 0). In the case of 
single-level metalization, from (39) and (49), 

= - (Ca + D)- 1 1 F n (x)\ n E n - f i *»(x) ^^ , (53) 

where 

_ M' n (0) . ,~as 

En - wW L (54) 

It should be noted from (33), (35), (37), and (38) that, as n->«>, 

En ~ 2 ' ,M„(0) ~ 1 + v ' { ^ 

In the case of double-level metalization, from (39) and (50), 

^2 {x> o) - ?p (x, 0) = *^-° - (Ca + D) - t F n (x)X n H n 
dy dy 2/i_i n=i 

- i o n (x) . , x ; , - £ *„(*) -%^ . ( 56 ) 

n =i sinh X„Li n =i i?M„(0) 

where 

*--nSr **»•»->■ (57) 

It again follows from (33) and (37) that, asn-*», 

tf»«-2. (58) 

For a periodic function /(x), of period L, we denote by c„[/] and 
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Fig. 4 — The splines used to represent the charge density on the electrodes. 
s„[_f] the cosine and sine Fourier coefficients of /: 

CnUl = \ /* f(x) cos (X n x)dx, s„[/] = | j* f{x) sin (\ n x)dx. (59) 

Then, from (52), the Fourier series for p(x) is 

1 I N 

p(s) = 2 C °M + E ^' c °tP>] 

+ E Cn[p ] + E PiCn[Pj]l COSX n X 

+ \sn[_p ~] + E P/9s[py]| sinX„x. (60) 
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In Appendix A we give c„[p y ] and s„[p>] for the various functions 
Pj(x). If we now equate (53) with (60), we can obtain the a„ and b n 
as linear functions of the />, in the case of single-level metalization: 

ao _ 2JEC-_AD) _ A j co[po] + £ p>CoM j f (61) 

a n =- itnL n (hi)l/lvKE n M n (0)l 

- | Cn[p ] + .E PiC,M } / (XA), (62) 

fon =- [^(Wl/E^^M^O)] 

- j «.[> J + .£ Pi*»CPi] 1 / (Xn-B n ). (63) 

Similarly, if (56) is equated with (60), we obtain the a n and b n in the 
case of double-level metalization : 

_ 2h- i(BC - AD) - apA _ ^fe-i 
a ° " (/i_0C - A (ft-i)C - A 



CoW+E PiCoCPi]}, (64) 



Q!n _ tnLnjhl) 

an ~ H n sinh (/i_i\„) i)\„i?»M„(0) 



&. =- 



( N 

~ UnlPal + E ViCnlvi] 
Bn £nLn(hl) 



(\ n H n ), (65) 



H n sinh (fc-iXn) vKH n M n (0) 

- J«J>,] + E P>Sn[P,][ /(Xn^n). (66) 

Equations (39) to (42) and either (49) or (50), with the a„ and 6 n 
denned by eqs. (61) to (63) or (64) to (66), respectively, define solutions 
to eqs. (21) to (25) which satisfy boundary conditions (11) or (12), 
(14) to (16), (26), and (27). They do not, however, assume the correct 
values on the electrodes at y = 0; i.e., (13) is not satisfied. These 
solutions depend linearly on the N unknown constants p j} I £ j £ N, 
and of course on the choice of functions p s (x) used to describe the 
charge density on the plates. Having picked the Pj(x) described 
earlier, Mj points are chosen on the jth electrode, x^, 1 ^ i ^ M jt with 

M = E M ,- £ N. 

y-i 
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Then the expression 

t Z LVo(.rP, 0) - Vtf (67) 



is minimized with respect to the p a . 

It has been observed that near the edge of an electrode, where the 
potential has square-root behavior, the fitting points x\ n should be 
spaced quadratically closer together as the edge of the electrode is 
approached. If x = e is an edge of the jth. plate, then we distributed 
the points near this edge by picking b t± e and an integer m < Mj/2 
and setting 

Zj ■ = b + (e - b) cos fc^j , (1 £ J £ m). (68) 

The points «,- were then used as the fitting points x* near the edge of 
the plate. Away from the edges of the plate, the fitting points were 
uniformly distributed. 

We have assumed the existence of a bounded solution \p(x, y) to the 
linearized problem. In Appendix B we show that if \p{x, y) is the true 
solution of the linearized problem and \f/ a (x, y) is one of our approximate 
solutions, then the error >f/(x, y) — \f/ a (x, y) is bounded at every point 
by the maximum error on the electrodes. Since the true solution is 
known on the electrodes, this provides us with a posteriori error bounds. 
We will make use of this important point later in evaluating the 
quality of our approximate solutions. 

The technique described in this section can be formulated in a 
rather general setting and, we believe, can be applied to many problems 
of interest in physics and engineering. It has been used by Morrison 
et al. in a study of microwave scattering by deformed raindrops. 11 
Assume that a problem can be separated into two parts: Input data 
and a governing system of partial differential equations (PDE's), with 
possible interface conditions, which determine the solution when given 
the input data. Further, assume that linear families of particular solu- 
tions to the PDE's can be found. For example, these may be con- 
structed by separation of variables, Fourier series, Green's Theorems, 
etc. Finally, assume that by linearly parameterizing some unknowns 
of the solution (for our problem, the charge distribution on the plates) 
we can obtain particular solutions to both the PDE's and the interface 
conditions. Then one could use some fitting procedure, a discrete 
least-squares fit, for example, to force the linear family of particular 
solutions to the governing system to have approximately the same 
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input data as the desired solution. This generates a solution of exactly 
the same governing system, but with input data which differs by a 
known, and hopefully small, amount from the desired input data. For 
all practical purposes, this gives an effective bound on the error in the 
computed approximate solution. For example, if the desired input is 
only known to 1 percent, because of experimental error in measuring 
it, then any solution generated by the above procedure which corre- 
sponds to the desired input data perturbed by at most 1 percent 
cannot, on the basis of comparing inputs, be distinguished from the 
true solution of the problem. 

Also, in many cases, one can use the Maximum Principle, conserva- 
tion of energy, or some other basic principle to give sharp, rigorous 
bounds on the error in such an approximate solution in terms of how 
well it satisfies the given input data. We do this for our problem in 
Appendix B. This is a very great improvement over the standard 
discretization methods for solving such problems. Those methods 
generally give an approximate solution to an approximate system of 
equations, but with exactly the given input data, with the result that 
it is very difficult to estimate reliably what the true error is for a given 
approximate solution. 

IV. THE POTENTIALS AND FIELDS IN SOME SPECIFIC CCD's 

Using the method described in Section III, we have evaluated approx- 
imately the solutions of eqs. (21) to (25) for a number of different plate 
configurations and design parameters, and we present some of these 
results graphically in this section. 

We have assumed in each case that the n-type substrate doping is 
N* D = 10 14 cm- 3 , that e 2 /e = 12, where e„ is the permittivity of free 
space, that ei/e 2 = J, and that Q(x) s 0, i.e., there is no trapped or 
implanted charge at the oxide-semiconductor interface. Then at 
T = 300°K, the Debye length is X = 0.415 ^m. In addition, we have 
used the factor (kT/e) = 0.025 V to convert dimensionless potentials 
to volts, and the factor (kT/e\ D ) = 600 V/cm to convert dimensionless 
fields to volts per centimeter. In each example involving a buried 
channel CCD, we assume that the acceptor number density in the 
p-layer is given by (1), [(6)] with (?, = 4.6 X 10 15 cm" 3 [C. = 46]. In 
each such case, this corresponds to an average number density of 
acceptor atoms of 2 X 10 15 cm -3 [see eq. (2.5) of I]. 

In I we investigated the effects of changing the p-layer doping and 
thickness, and so here we concentrate mainly on the effects of gap 
width, plate potential, and the separation of the levels of metalization. 
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Fig. 5 — The channel potential <p* plotted as a function of x* for a three-phase 
buried channel CCD. The 45-/im-wide electrodes are at 0, —5, and —10 V; the gaps 
are 5 Mm wide; /if = 0.1 urn; h% - hf = 5 nm; and CT = 4.6 X 10 18 cm -3 . The 
dashed curve is for no surface charge implanted in the gaps ; the solid curve is for 
pi/e = 0.8 X 10 12 cm" 2 implanted in the gaps. 



In Figs. 5 and 6 we show some properties of a three-phase buried 
channel CCD with single-level metalization. The electrodes are 45 /*m 
wide, and the gaps between them are 5 /urn wide. The p-layer is 
5 nm thick (h* 2 — h\ = 5 /xm), and the oxide layer is 0.1 pm thick 
(hi =0.1 jum). The region y < is assumed to be filled with SiO 2. The 
potentials on the electrodes are 0, —5, and —10 V, as shown. The 
dashed curve in Fig. 5 shows the channel potential <p* (that is, the 
value of the potential at the potential minimum in the p-layer) as a 
function of x* when there is no implanted surface charge in the gaps 
between electrodes \j>g(x*) m 0]. This curve illustrates one early 
difficulty encountered in the design of buried channel CCD's, namely 
the large potential well under the gap between the plates. A CCD 
with almost these same parameters was constructed 5 and did not work 
because of the variable amounts of charge trapped in these wells. In 
the remainder of this section we discuss a number of possible methods 
of eliminating this potential well in the gaps between the plates. 
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Fig. 6— The channel field -d<p*/dx* plotted as a function of x* for the CCD of 
Fig. 5 with pj/e = 0.8 X 10 12 cm -2 implanted in the gaps. 



An operating buried channel CCD has been reported 12 in which the 
gaps between the electrodes have been filled with a resistive material 
so that the potential drop between the electrodes is essentially linear. 
This CCD was also discussed in I, and it was shown there that the 
potential wells are eliminated. Another technique for eliminating the 
potential wells is to implant a layer of positive surface charge in the 
gap between the electrodes. (Other schemes for eliminating this 
problem are discussed in the literature. 13 ) The solid curve in Fig. 5 
shows the channel potential in the same three-phase CCD after a uni- 
form surface charge density, p* g (x*)/e = 0.8 X 10 12 cm" 2 \_p (x) = 578], 
has been implanted in the gaps between the electrodes. Note that 



r h ' 
Jh\ 



N\{y*)dy* = 10 12 cnr 



This technique should also eliminate the potential wells under the 
gaps. In Fig. 6 we plot the channel field E* x = - d<p*/dx* (that is, the 
field at the potential minimum in the p-layer) as a function of x*. 
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Fig. 7 — The channel potential in a buried channel CCD with double-level metaliza- 
tion. The lower level electrodes are 10 nm wide with 5-/im gaps and are at potentials 
of —5 and —7 V. The upper level is a single electrode at —5 V; /if = 0.1 ^m; 
/i* — hf = 5 nm ; C? = 4.6 X 10 16 cm -3 ; and the separation of the metalization levels 
is h*_i = 0.1, 0.5, and 1.0 ftm. 

This shows that there are substantial fields in the gap between the — 5 
and —7 V electrodes, but the field penetration under the electrodes is 
not too good. 

In Figs. 7 to 9 we investigate the possibility of eliminating the poten- 
tial wells by the use of double-level metalization, in which the upper 
level of metal is a single continuous piece covering the entire channel 
and has a dc potential applied to it. The presence or absence of the 
potential wells is mainly a local phenomenon and can be studied by 
considering just two adjacent electrodes in the lower level of metaliza- 
tion. Thus, in the interests of economy we consider a model CCD in 
which alternate electrodes of the lower level are at the same voltage. 
These plates are 10 m^i wide and the gaps between them are 5 /xm 
wide. The oxide layer between the first level of electrodes and the 
p-layer is 0.1 /an thick (h\ =0.1 fim) and the p-layer is 5 /jm thick 
(h* 2 — K = 5 /zm). For the CCD of Fig. 7, the electrodes on the lower 



1784 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1973 

I -5V 



-22.50 



-23.75 - 



-25.00 - 



-27.50 - 



-28.75 - 



30.00 - 



-31.25 - 



-32.50 - 



-33.75 




12 14 16 



28 30 



Fig. 8— The channel potential in the CCD of Fig. 7 with the separation of the 
metalization levels held fixed at h^ = 0.5 urn and the potential of the center, lower- 
level electrode taking the values —7, —10, —15, and -20 V. 



level are at a potential of -5 and -7 V, and the upper level consists 
of a single electrode, covering the whole device, which is at a potential 
of — 5 V. In Fig. 7 we plot the channel potential for three different 
separations of the levels of metalization, /i*_i = 0.1, 0.5, and 1 nm. Even 
with h'-! =0.1 /xm, there is still a very slight well in the gaps. 

We next see what happens if we hold the separation of the levels of 
metalization at h*- x = 0.5 m and change the voltage, v , on the middle 
electrode in the lower level. In Fig. 8 we plot the resulting channel 
potential, <p*, as a function of x* for v = —7, -10, -15, and —20 V. 
From these graphs we see that a potential difference of 15 V between 
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Fig. 9 — The channel potential in the CCD of Fig. 7 with the separation of the 
metalization levels held fixed at /i*i =0.1 /im and the potential of the center, lower- 
level electrode taking the values — 7 and — 10 V. 



neighboring electrodes on the lower level will insure the absence of 
potential wells in the gaps. 

In Fig. 9 we plot the channel potential for the same device but 
with /i*_! =0.1 nm and for v = — 7 and —10 V. For this small separa- 
tion of the levels of metalization, a 5-V potential difference between 
neighboring electrodes eliminates potential wells in the gaps. 

In Figs. 10 and 11 we show some effects of gap width. We consider 
first a buried channel CCD with double-level metalization. The p-layer 
is 5 ^m thick (h\ — h\ = 5 //m) , the oxide layer between the first level 
of electrodes and the p-layer is 0.1 ^m thick {h\ =0.1 /xm), and the 
layer between the two levels of electrodes is 0.5 nm (h*-i = 0.5 Mm). 
The upper level of electrodes consists of a single electrode at a potential 
of —5 V. The lower level consists of electrodes 10 fiin wide and, as 
shown in Fig. 10, they are alternatively at potentials of —5 and —7 V. 
Curves of the channel potentials are plotted for three different gap 
widths between plates: 5, 1, and nm. (The 0-gap curve was calculated 
by the methods of I.) The x* scale for the three curves are different, but 
are chosen so the centers of the gaps coincide. With 5-/xm-wide gaps, 
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Fig. 10— The channel potential in the CCD of Fig. 7 with the separation of the 
metalization levels held fixed at k*, =0.5 /xm, and the gap width taking the values 
5, 1, and /xm. 

there are large potential wells, as we saw in Figs. 7 and 8. However, by 
reducing the gap width to 1 fim, the potential wells are essentially 
eliminated. The curve for zero electrode separation is included to show 
that it is a good approximation to the channel potential in cases of 
small electrode separation. 

Finally, in Fig. 11, we plot the potential along the oxide-semicon- 
ductor interface (y* = h\) for two surface CCD's. In each case, the 
oxide layer is 0.1 jum thick (h\ = 0.1 /im) and the region y* < is 
assumed to be filled with Si0 2 . Also, in both cases, the electrodes are 
10 mhi wide and are held at alternate potentials of -5 and -7 V. In 
one case the gap between electrodes is 1 *im, while in the other case 
there are no gaps between electrodes. (The zero gap curve was calcu- 
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Fig. 11 — The potential f*{x*, /if) plotted as a function of x* for a surface CCD 
with single-level metalization ; h* = 0.1 ^m; 10-pm-wide electrodes held at —5 and 
—7 V; and with gaps between electrodes of 1 and tun. 

lated by the method of I.) Again the x* scales differ, but the centers of 
the gaps coincide. Except in the region between the plates, the two 
curves coincide closely. 



V. COMMENTS ON ACCURACY AND COST 

We considered in some detail in I how well the solutions of the 
linearized equations (21) to (25) approximate the solutions of the 
nonlinear equations (7) to (10). It was shown there and in Ref. 4 that 

as long as max fa(x, hi) ^ — 160, and I fa(x, h 3 ) + 1 1 < 10, ^ x ^ L, 

osxgz- 

[|^a(z, h 2 ) + 1| < 10 for surface devices] then the solution of the 
linearized problem approximates the solution of the nonlinear problem 
to within several percent in the p-layer for buried channel devices 
(near the oxide-semiconductor interface for surface devices). In addi- 
tion, we are interested in examining the accuracy of the approximate 
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solutions of the linearized equations. As we stated in Section III (and 
prove in Appendix B), the difference between the exact solution and 
the approximate solution, 8*(x*, y*) = **(**, y*) - \p a *(x*, y*), is 
bounded everywhere by its maximum value on the electrodes. 

As typical examples, consider the curves in Figs. 10 and 11. For 
all three curves in Fig. 10, we found that fa(x, fa) ^—345 and 
\faix, fa) + 1 1 < 2.75 for ^ x g, L. For the two curves of Fig. 11, 
we found that if/i(x, fa) < - 186 and | fa(x, fa) + 1 1 ^ 1 for ^ x ^ L. 
Furthermore, for the curves of Fig. 10, a search of the electrodes showed 
that, for 5-Mm gaps between electrodes, max 1 8*(x*, y*) | ^ 0.29 V 
and, for l-/zm gaps, max | 5*(x*, y*) \ ^0.18 V. (The no-gaps curve was 
calculated by the methods of I.) These correspond to maximum errors 
of the channel potential of 1.3 and 0.83 percent respectively. For the 
curve for the surface CCD with l-/xm gaps between the electrodes a 
search of the electrodes showed that max 1 8*(x*, y*) | ^ 0.4 V. (The 
no-gaps curve was again calculated by the methods of I.) This corre- 
sponds to a maximum percentage error of 0.86 percent. 

For the curves of the buried channel CCD of Fig. 10, we have 
undoubtedly overestimated the error in the channel for the following 
reasons. It can be shown that the error in each coefficient a n and b n 
appearing in (34) can be expressed as an integral over the electrodes of 
the error 8(x, y) times a weight function. The sign of 8(x, y) oscillates 
on the electrodes, and so one would expect the error in the lower-order 
coefficients to be quite small. Furthermore, an examination of (34) to 
(38) and (40) shows that, for the parameters involved, only the first 
ten terms in (40) contribute significantly to the channel potential. 

To present some idea of the cost of running these programs, the 
calculation of the solution for the case of 5-//m gaps in the curves of 
Fig. 10 took 253 seconds and used 40 K of core, the case of l-/im gaps 
took 258 seconds and 40 K of core. Calculation of the solution for the 
l-,um gap case of Fig. 11 took 263 seconds and used 40 K of core. By 
comparison, the two corresponding no-gap solutions, obtained by the 
methods of I, took 16 seconds and 32 seconds, respectively, and both 
solutions required 40 K of core. 
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APPENDIX A 



This appendix lists the coefficients of the Fourier series of the 
various splines used in approximating the charge density on the elec- 
trodes. We begin by defining the function 



u(x, h) = 



V n ' (69) 



0, h ^ x ^ L. 



Outside the interval S % ^ L, this function is defined by periodicity, 
u(x, h) = u(x + L, h). The Fourier coefficients of u(x, h) are 

2 

a n (h) = c n \_u] = j-rr^ (1 — cos \ n h), 

" 2 (70 > 

b n (h) = s n [w] = Yj—z (\ n h — sin X n h), 

where the notation c„[u] and s„[w] is defined in (58), X„ = 2mr/L, 
and n = 0, 1, 2, • • ■. Note that 

a (h) = ^ , b (h) = 0. (71) 

The triangular splines can all be expressed in terms of u(x, h), and 
their Fourier coefficients are simple linear functions of the coefficients 
a n (h) and b n (h) defined in eq. (70). Thus (see Fig. 4a): 

p (x, h) = u(x, h) + u(L - x, h), (72) 

and for n = 0, 1 , 2, • • • , 

c„[p ] = 2a n (h), s„[p ] = 0. (73) 

Similarly (see Figs. 4b and 4c), we have the end splines 

pt(x; x , h) = u(x — x , h), p r (x; xo, h) = u(x — x, h), (74) 

and for n = 0, 1, 2, • • • 

c*\j>i] = (cos \ n x )a n (h) — (sin \ n x )b n (h), 
s n£pi~] = (sin \ n xo)a„(h) + (cos X„xo)b n (h), 



and 



c «[Pr] = (cos \ n x )a n (h) + (sin \ n x )b n (h), 
s„[p r ] = (sin \„x )a n (h) — (cos \ n Xo)b n (h). 



(75) 
(76) 
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Finally (see Fig. 4f), 

t(x; x , h h h 2 ) = u(x — x, hi) + u(x — x 0) h 2 ), x ^ a,- (mod L) 

= 1, x = £ (niod L) (77) 

and 

c n [t] = cos X„a:o{a n (^x) + a„(h 2 )} + sin \ n x {b n (hi) - b n (h 2 )}, 
s n lQ = sin\„a;o{a„(/ii) + a n (h 2 )\ — cos \„x {b n (hi) — b n (h 2 )}. 

The coefficients of the Fourier series of the singular edge splines were 
calculated as follows. We define (see Figs. 4d and 4e) 

f = (x — x )~ h — h~l, x < x ^ x + h, 
s ( (x; x 0t h) = \ (79) 

[ = 0, ^ x ^ x , x + h < x ^ L, 

and 

s T (x; Xo, h) = st(2x ~ x;x ,h). (80) 

Then after an integration by parts 

c n |>] = 

4 f ,. . , . , N sin X„(z + h) — sin \„x 
= z |W cos X n (x + h) 2X^P 

/xo+h "J 

(x — Xo)* sin \ n xdx , (81) 

s*M = 

4 f, . . . , . ,v . cos X„(a;o + h) — cos \ n x 
= Z [ h Sm (X ° 2XTP 

/xo+h "I 

(z — Zo)* cos X n zd:r. . 

The integrals on the right of (81) were evaluated by quadratures using 
Filon's method. 14 Similarly, 

Cn[Sr] = 

4 T, , . , , N sin X„x — sin \„(x — h) 
= z \h* cos X n (x -h) 2JOi 

— X„ / (x — x)> sin \ n xdx , (82) 

Jx a -h J 

Sn[S r 2 = 

4 [".. . . , IN . cos X„a:o — cos X„(x — h) 
= j \h* sin X n (x — h) H 2x~p 



+ X„ / (zo — z) J cos X„zd:r . 
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APPENDIX B 

We assume that a solution of equations (21) to (25) satisfying 
boundary and interface conditions (11) or (12) and (13) to (16), (18), 
(26), and (27) exists and is bounded. As before, we denote this solution 
by $( x i y)> we denote our approximate solution by ^ a (x, y), and we 
define 

ffo y) = f(x, y) - p(x, y). (83) 

By \p a , \pa, and £ a , (0 ^ a ^ 4), we mean \J/, \p a , and £ restricted to the 
various subdomains. 

From their construction, the approximate solutions satisfy the same 
equations and boundary and interface conditions as the exact solution, 
except they do not assume the correct values on the electrodes. 
Consequently 

v««.(*, y) = o, (o ^ a ^ 3), v*Ux, v) - Ux, y) ="0. (84) 

Also £a(x, y) satisfies either boundary condition (11) in the case of 
single-level metalization or 

b(x, -/>_,) = (85) 

in the case of double-level metalization. In addition, £o(x, 0) = £ i(x, 0), 
^ x ^ L, and $ a (s, y) satisfies (14) with Pg = 0, (15) with Q(x) = 0, 
(16), (18), (26), and (27). 

We now outline a proof that if M = sup j £o(z, 0) | , where E denotes 

the electrodes, then | £(x, y)\ ^ M for all (x, y). The plan of the proof 
is to show that £o(x, y) is bounded both above and below by its maxi- 
mum and minimum on (0 ^ x ^ L, y = 0); that £ a (x, y) is bounded 
above and below by its maximum and minimum on either (0 ^ x ^ L, 
y = /i a _i) or (0 ^ x ^ L, y = h a ), (a = 1, 2, 3), where we define 
ho = 0; and that ^(x, y) is bounded above and below by its maximum 
and minimum on (0 ^ x ^ L, y = h 3 ) . Then we show that the global 
maximum and minimum must occur on the electrodes. 

First consider Ho(x, y) in the case of single-level metalization. Then 
£o(x, y) is harmonic in the strip S = (0 ^ x ^ L, — °o < y ^ 0), and, 
by the Phragmen-Lindelof theorem [Ref. 15, corollary to theorem 19, 
Chapter 2, with w(x, y, = 1 — y], £o(x, y) is bounded in S , both above 
and below, by its values on the lines (x = 0, — « < y ^ 0), (y ^ x ^ L, 
y = 0), and (x = L, - co < y ^ 0). Let m = inf f (z, 0) and 

Mo = sup £o(x, 0). [Note that since £(.r, 0) is continuous, there exist 

points ^ x m , x M ^ L such that m = £o(x m , 0), M = £o(x M , 0).] 



1792 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1973 

Further, we must have 

Zo(x, y) = h@o + L Q3„ cos \ n x 4- 7» sin \ n x)e K ^, (86) 



and hence 



1 1 f L 

lim £o(0,2/) = lim b(L,y) =^Pq= t / £0(2, 0)(fe. 



(87) 



From (87) we can conclude that 

mo ^ lim £ (0, y) = lim £ (L, y) ^ M . (88) 

Now if Mo is not the maximum value of £o(z, ?/), from what we have 
just shown, and from the periodicity of £(x, y) in x, this maximum 
value must be assumed at two points (0, y ) and (L, y ), with 

— °° < ?/o < 0. Further, the outward directed normal derivative at 
these points must be positive (Ref. 15, theorem 8, Chapter 2); that is, 

— (d(n/d£)(0, yo) > 0, (dto/dxXL, y ) > 0. However, from periodicity, 
(db/dx)(Q, j/o) - (d£a/dx)(L, y ), which is a contradiction. Hence M 
is the maximum value of £o(x, y) in S . The same reasoning applied to 

— £ (x, y) shows that m is the minimum value in S . 

In the case of double-level metalization, the maximum principle for 
harmonic functions (Ref. 15, theorem 2, Chapter 2), plus the boundary 
condition fa(x, —h-i) = 0, implies that £ (&, y) is bounded everywhere 
in (0 ^ x ^ L) X ( — ft_i ^ 2/ ^ 0), both above and below, by its 
values on the sides (x = 0, -/i_i ^ y ^ 0), (0 £ 3 £ L, 0), and 
(x = L, —h-i^y^0). Then the same reasoning as in single-level 
metalization shows that £ (x, y) achieves its maximum and minimum 
on (0 ^ x ^ L, y = 0). 

Essentially the same arguments used in the double-level metalization 
case can be used to show that £„(x, y) (a = 1, 2, 3) must achieve both 
its maximum and minimum either on the line (0 ^ x ^ L, y = /i„_i) 
or (0 £ x £ L, y = h a ), where we define h Q = 0. 

The Phragmen-Lindelof theorem can be applied to £ 4 (x, y) in 
S 4 = (0 ^ x ^ L) X (h 3 ^ y < <*) to show that £ 4 (x, y) is bounded 
everywhere in «S 4 both above and below by its values on (x = 0, 
h 3 ^ y < * ), (0 ^ x ^ L, ?/ = /i 3 ), and (x = L, h 3 ^ y < «). Making 
use of the boundary condition f 4 («, «> ) = 0, the same arguments used 
in the case of single-level metalization for £ (x, ?/) show that £ 4 (x, y) 
achieves its maximum and minimum values on (0 & X ^ L, y ™ »$). 

We have shown that £(z, y) must assume its maximum and minimum 
values at points on the lines (0 S x S L, y = h a ), ^ a ^ 3. Suppose 
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that £(x, y) is a global maximum at the point P = (z, h a ). Then clearly 
\p a (x, y) and \p a +i(x, y) take on their maximum values at P. Conse- 
quently, their exterior normal derivatives at P must be positive (Ref. 
15, theorem 8, Chapter 2), i.e., (d+Jdy)(P) > 0, - (#«+i/dy)(P) > 0. 
However, if P is not a point on an electrode, it follows from inter- 
face conditions (14), (15) with Q(x) = 0, (16), or (26) that either 
n(ty a /dy)(P) = (d+ a+ i/dy(P) or (ty a /dy)(P) = (d+ a+1 /dy)(P), 
which is a contradiction. Consequently £(z, y) achieves its maximum 
Mo on an electrode. The same argument applied to — £(x, y) shows 
that £(z, y) also achieves its minimum m on an electrode. If we set 
M = max(|m |, | il/ |), then we have shown that \%(z, y)\ ^ M. 
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